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Abstract 


The water level reduction in Lake Urmia and its effects on the surrounding environment has been 
among the important national and international challenges in the past two decades. Therefore, this 
study aimed to investigate the current status and predict the future state of land use in the Gadarchai 
watershed, located in West Azerbaijan Province, which is one of the important sub-watersheds of 
the Lake Urmia basin. For this purpose, Sentinel-2 satellite images for 2016, 2020, and 2022 were 
obtained from the European Union's Copernicus website. Then, preprocessing methods were applied 
in various software environments, and the relevant images were sent to the eCognition software 
environment. In this environment, using various basic and object-oriented methods (especially 
segmentation and production of different coefficient layers), the nearest neighbor classification 
method was implemented, and land use maps were produced. Finally, the Markov chain model was 
used to predict changes in land use in future years. To verify the accuracy of the Markov chain 
model, the predicted land use change map for 2022 was compared with the 2022 classification map. 
The research results showed that with the application of basic methods, especially nearest neighbor 
classification, it is possible to produce land use maps with high accuracy (90% kappa coefficient). 
Also, by applying the Markov model, land use change maps with an acceptable accuracy level 
(around 80%) are possible. The final results indicate that by the year 2028, agricultural land use 
(13.89%), dry farming (14.1%), residential areas (0.33%), and salt pans of Lake Urmia (26%) will 
increase. It should be noted that the soil use class will decrease by 10.26%, and pastures will decrease 
by 5.3596. Overall, the final models demonstrate the high accuracy of basic and object-oriented 
methods and the suitable performance of the Markov model in the process of studying land use 
changes. 
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حکیده 


کاهش سطح تراز آبی دریاچه ارومیه و اثرات آن بر محیط پیرامون دریاچه از موضوعات و چالش‌های مهم ملی و بین‌المللی 
در دو دهه اخیر بوده است. بر اساس مطالعات صورت گرفته یکی از مهم‌ترین عامل اثرگذار بر این روند تغییرات اراضی. 
به‌ویژه کشاورزی بوده است. بر همین اساس هدف تحقیق حاضر بررسی وضع فعلی و پیش‌بینی وضعیت آتی کاربری 
اراضی حوزه آبخیز گدارچای واقع در استان آذربایجان غربی یکی از زیر حوضه‌های مهم حوضه آبریز دریاچه ارومیه است. 
به همین منظور ابتدا تصاویر ماهواره‌ای سنتینل-۲ برای سال‌های میلادی ۲۰۱۳۰۲۰۲۰ و ۲۰۲۲ از سایت کپورنیک اتحادیه 
اروپا دریافت گردید. سپس روش‌های پیش‌پردازش در محیط نرم‌افزارهای مختلف اعمال و تصاویر مربوطه به محیط 
نرم‌افزار 0 ارسال شد. در این محیط با استفاده از روش‌های مختلف دانش ub‏ و شی گرا (به‌ویژه سگمنت 
سازی و تولید لایه‌های ضرایب مختلف) روش طبقه‌بندی نزدیک‌ترین همسایگی اجرا و نقشه‌های کاربری اراضی تولید شد. 
درنهایت از مدل زنجیره مارکوف به‌منظور پیش‌بینی تغیبرات کاربری اراضی در سال‌های آتی استفاده گردید. برای بررسی 
صحت مدل CA‏ مارکوف. ابتدا نقشه تغییرات پیش‌بینی شده سال ۲۰۲۲ با نقشه طبقه‌بندی ۲۰۲۲ صحت‌سنجی شد. quU‏ 
تحقیق نشان داد که با کاربرد روش‌های دانش ah‏ به‌ویژه طبقه‌بندی نزدیک‌ترین همسایگی امکان تولید نقشه‌های 
کاربری اراضی با دقت بالا (ضریب ٩۰ LIS‏ درصد) امکان‌پذیر است. ضمناً با اعمال مدل مارکوف نقشه‌های تغییرات 
کاربری اراضی با دقت قابل‌قبول (در حد ۸۰ درصد) امکان‌پذیر است. نتایج نهایی مبین این واقعیت است که تا سال ۲۰۲۸ 
میلادی کاربری اراضی کشاورزی (۱۳۸۹) کشت دیم (M‏ محدوده‌های مسکونی (۰/۳۳) و عرصه‌های Sai‏ دریاچه 
ارومیه حدود (VV)‏ درصد افزایش خواهد داشت. یادآور می گردد کلاس کاربری خاک در حد (۱۰/۲7) و مراتع به میزان 
(0/Y0)‏ درصد کاهش خواهند داشت. درمجموع مدل‌های نهایی مبین دقت بالایی روش‌های دانش‌پایه و شیءگرا و کارایی 
مناسب مدل مارکوف در روند مطالعه تغییرات کاربری اراضی هستند. 
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مقد مه 

یکی از مهم‌ترین عامل‌ها برای درک تغییر محیطی در تمام مقیاس‌های زمانی-مکانی تغییر کاربری‌زمین 
است )2015 (Kumar, Bhaskar, & Padmakumari,‏ تغییر کاربری زمین با عامل‌هایی مانند افزایش ارتفاع 
در حوزه‌های آبخیز افزايش ظرفیت سیلابی رودخانه‌ها و افزايش ارتفاع آب بر حجم رواناب در 
رودخانه تأثیر خواهد داشت )2017 (Moe, et aL,‏ آگاهی از تغییر کاربری زمین یکی از مهم‌ترین 
پارامترهای تصمیم‌گیری برای برنامه‌ریزی زمین در مقیاس جهانی و ناحیه‌ای است. پیشرفت در درک 
حالت‌های ممکن تغییر کاربری‌زمین می‌تواند در شناسایی و تشخیص فرآیندهای زمین (سیلاب و...) در آینده 
مناسب باشد )2019 (Armenteras, Murcia, González, Barón, & Arias,‏ بررسی تغییرات می تواند 
به درک فرآیندهای انسانی کمک کرده تا بتوان از میزان تخریب ایجاد شده جلوگیری شود (Donnay,‏ 
Longley, & Barnsley, 2001)‏ تهیه نقشه تغییرات کاربری اراضی. بخش زیادی از اطلاعات موردنیاز 
متخصصان و برنامه‌ریزان توسعه شهری و روستایی را فراهم می‌کند. تهیه این اطلاعات در سطح وسیع و با 
دقت بالا از طریق انجام عملیات زمینی امری هزینه‌بر زمان‌بر و در مواردی غیرممکن است. امروزه فناوری 
سنجش از راه دور راهکاری ارزشمند در جهت شناسایی. تشخیص, پایش. و پهنه‌بندی منابع‌طبیعی به‌ویژه در 
مورد روند ag‏ نقشه‌های کاربری اراضی» در مناطق مختلف جهان به‌صورت علمی مورد استناد قرار Jal‏ 
گرفت (رسولی»۱۳۷۸). 

منطقه‌های وسیع حوزه‌های آبخیز توسط داده‌های سنجش‌ازدور پوشش داده می‌شوند و همچنین پوشش 
زمانی زیاد این داده‌ها اطلاعات باارزشی از مکان, فرآیند نسبت» روند طبیعی و الگوی تغییر پوشش کاربری 
زمین فراهم خواهد کرد )2019 (Munthali, Botai, Davis, Ade La Abiodun,‏ برای شناخت TIG‏ 
کاربری‌زمین و مکانیزم‌های پیچیده مدل‌ها و تحلیل حالت‌های ممکن ابزاری توانا است. برخی از این 
مدل‌های شبیه‌سازی ممکن است به تصمیم‌گیری در مورد پیش‌بینی حالت‌های که در آینده اتفاق می‌افتد 
کمک کند. ازجمله مدل‌های که می‌توان نام oy‏ مدل‌ها اتوماتای سلولی و زنجیره مارکوف است 
(Samie, Deng, Jia, & Chen, 2017)‏ مدل‌سازی مکانی تغییر کاربری‌زمین به اطلاعاتی در مورد اندازه 
تغییر و این‌که تغییر در کجا اتفاق خواهد افتاد (مقدار و مکان) نیازمند می‌باشد. مدل زنجیره مارکوف 
اطلاعاتی از مکان نمی‌دهد و تنها پویایی زمانی را به‌دست می‌آورد اما در مقایسه با مدل زنجیره 
ab Syl‏ مدل اتوماتای سلولی CA!‏ دارای مولفه مکانی می‌باشد» که با قانون‌های خاصی از سلول‌های 
همسایه تغییر در آینده به‌دست می‌آورد. خروجی مدل مارکوف از نظر ماهیت. مکانی نیست؛ یعنی در آن از 
موقعیت جغرافیایی کاربری‌های زمین هیچ دانش و آگاهی وجود ندارد. زمان 1+1 برای پیش‌بینی موقعیت 
کاربری‌ها. به همراه مدل مارکوف روش سلول‌های خودکار به کار گرفته می‌شود. مدل CA‏ مارکوف 


1. Cellular Automata 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف ۲۱۹ 


دانش کاربر و مولفه‌های مجاورت مکانی را به توزیع مکانی به مدل زنجیره مارکوف JU al‏ تبدیل 
مفید است چون روند تغییر کاربری زمین و حالت ممکن در نظر می‌گیرد (Mujiono, Harmantyo,‏ 
Rukmana, & Nadia, 2017)‏ به‌طورکلی توانایی شبیه‌سازی و مدل‌سازی حالت‌های ممکن مختلف با 
به کاربردن روش‌های نوين در پردازش تصویرهای Slo yl gale‏ 3 مدل‌های کار dal‏ در تحلیل‌همای مکانی. 
امکان پیش‌بینی و بررسی تغییر زمانی - مکانی کاربری زمین را می‌دهد & (Ebrahimy, Rasuly,‏ 
Abmadpour, 2019)‏ برای استخراج اطلاعات کاربردی طبقه‌بندی تصاویر ماهواره‌ای یکی از مهم‌ترین 
روش‌ها محسوب می‌شود )2010 (Blaschke,‏ در حال حاضر طبقه‌بندی تصاویر PX‏ ماهواره‌ای L‏ سه 
NS‏ روش جدید شی ء گرا برای تجزیه‌وتحلیل lam‏ استفاده می‌شود که قطعه‌بندی سطح onl‏ و بدون 
دانش را با روش‌های طبقه‌بندی فازی سطح YU‏ و مبتنی بر دانش ادغام می‌کند. این روش جدید از طریق 
نرم‌افزار Cognition‏ که شامل محیطی شىء گرا برای طبقه‌بندی تصاویر ماهواره‌ای بود» طراحی ido‏ 

(Tzotsos & Argialas, 2008)‏ در روش طبقه‌بندی شیءگرا اطلاعات ab‏ شکل و بافت‌ها را به‌طور 
هم‌زمان در نظر می گیرد. معمولا خطاهای متعددی در روش‌های طبقه‌بندی مبتنی بر پیکسل قابل‌انتظار است 
که با توجه به قدرت تفکیک مکانی سنجنده و ple‏ عوامل متأثرکننده باید شناسایی شوند؛ اما مجموعه‌ای از 
پیکسل‌های alts bud‏ عارضه مشیعصی در روش‌های شیءگرا OUS‏ می‌دهند؛ cul plo‏ مجموغه پیکسل‌های 
مشابه تحت یک عنوان شناسایی و سپس عملیات طبقه‌بندی و فرآیندهای تقسیم‌بندی اعمال می‌گردد 
(رسولی»۱۳۷۸). در این بخش از تحقیق به بررسی و مرور نتایج انجام شده» به gas AX‏ 45 بارز در خارج از 
کشور و ایران اشاره می‌گردد. )2019 L (Munthali, Botai, Davis, Ade La Abiodun,‏ استفاده از شبکه 
و سنجش‌ازدور تلاش Jis‏ تغییرات کاربری اراضی را برای بازه‌های زمان Yo‏ نا yeyo‏ بررسی کنند و از 
صحت Lis — 3 T‏ برای ارزیابی دقت نقشه‌ای Mg‏ شده استفاده کردند. نتایج نشان داد کاربری‌های 
مناطق مسکونی پهنه‌های آبی و زمین‌های pl‏ افزایش Cole‏ خواهند داشت و جنگل‌ها و زمین‌های 
کشاورزی کاهش مساحت خواهند داشت و همچنین بیشترین میزان تغییر مربوط به مناطق جنگلی به مقدار 
۶ درصد خواهد بود. همچنین در تحقیق دیگری )2021 *(Rasouli, Asgarova, & Safarov,‏ به بررسی 


تغییرات ۳۰ ساله LU/LC2‏ در منطقه قره‌باغ با استفاده از روش‌های تحلیل تصاویر ماهواره‌ای مبتنی برشی۳ 


1. Geographic Information System 
2. Land Use/ Land Cover 
3. Object-Based Image Analysis 
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(OBIA)‏ پرداختند. نتایج۱ (DTD‏ نشان داد که از سال ۲۰۱۹ تا ۲۰۲۱ تغییرات LC‏ کاملاً معنادار است. 
کاهش پوش شگیاهی AVY)‏ افزايش سطح بدون پوش شگیاهی (۱۱/۸/) و بیشترین تغییرات در اراضی 
تخریب شده با مساحت ۵/۱ کیلومترمربع YVA)‏ مشاهده شد. همچنین مساحت اراضی فرسوده ۳۵ درصد 
و اراضی بایر ۶/۶ درصد افزایش al,‏ است. دقت IS‏ ۰/۹۵ و LIS‏ ۰/۹۳ تغییرات قابل‌توجهی را در 
dab LC/LU‏ کرد. )2022( Igué, , Houndagba, Gaiser, & Stahr,‏ با استفاده از روش I Se oh‏ نقشه 
کاربری اراضی حوزه‌های اویم و اونین با دقت کلی ۸۷ درصد تهیه شد و به‌عنوان کارآمدترین روش معرفی 
گردید. همچنین در Ghafari, Moradi,& Modarres, (2018) olp!‏ « با استفاده از تصاویر TM2 cJ‏ 
سال ۱۹۸۵ ; OLI Y‏ سال ۲۰۱۵ به مقایسه روش‌های طبقه‌بندی شیءگرا و پیکسل پایه» در دشت‌های 
اصفهان- برخوار» نجف‌آباد و چادگان و بررسی تغییر کاربری اراضی در طول دوره آماری پرداخته‌اند. 
طبقه‌بندی کاربری اراضی شامل قطعه‌بندی داده‌های تصویری با استفاده نرم‌افزار eCognition‏ انجام گرفت. 
نتایج نماپانگر Kul‏ طبقه‌بتدی شیءگرا با صحت کی بالای ٩۰‏ درصد نسبت به طبقه‌بندی پیکسل‌پایه از 
دقت بالاتری برخوردار بوده است. )2020( Abiyat, Attar Roshan, & Abiyat,‏ در پیش‌بینی تغییرات 
بوشش گیاهی در ارتباط با تغییرات کازبری اراضی با استفاده از مدل 1۷4و زنجره CA‏ مار کرت در شهر 
اهواز از طبقه‌بندی شیءگرا استفاده کردند و نتایج نشان داد که پوششگیاهی روند کاهشی» و بیشترین 
تغییرات مربوط به تبدیل پوشش گیاهی به نواحی ساخته شده است و کمترین تغییرات مربوط به تبدیل 
پوشش گیاهی به پهنه‌آبی است و نتایج پیش‌بینی تغییرات تا سال ۲۰۲۹ نشان‌دهنده کاهش سطح پوشش ALS‏ 
می‌باشد. )2023( L Abedini, Pasban, & nezafat taklhe,‏ استفاده از روش شیءگرا نقشه کاربری اراضی 
حوزه آبخیز نیرچای در استان اردبیل تهیه کردند. با استفاده از طبقه‌بندی الگوریتم نزدیک‌ترین همسایه 
dS. uS‏ صحت طبقه‌بندی بررسی شد. نتایج OUS‏ داد که با استفاده از روش شیءگرا. صحت کلی ۹٩‏ 
درصد و ضریب UIS‏ ۰/۸۸ درصد در طبقه‌بندی به دست آمده است. علیرغم ارزشمند بودن تحقیقات انجام 
شده در مبحث بررسی تغییرات کاربری اراضی» در این تحقیق با دسترسی به داده‌های جدید فناوری 
سنجش زدور تصاویر Sentinel2‏ با SSH‏ مکانی ۱۰ متر و با اعمال روش‌های پردازش تصاویر ماهواره‌ای 
دانش‌پایه پرداخته است. از اهداف اصلی این تحقیق نظارت بر تغییرات زمانی-مکانی حوزه آبخیز گدارچای- 
eels‏ ارو کے کور Y EY‏ و ری کاربری cas‏ سال VA‏ توق dea d‏ فلز tidie‏ 
تکنیک پردازش تصویر سنجش‌ازدور با استفاده از تمام تصاویر موجود برای تشخیص تغییرات کاربری و 


پیش‌بینی آن برای آینده استفاده می‌شود؛ و درنهایت با استفاده از نتایج مدل» کاربری اراضی آینده پیش‌بینی 


1. Dynamic Thresholds Indexing 
2. Thematic Mapper 

3. Operational Land Imager 

4. Land Change Modeler 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف YYY‏ 


مدل‌های نهایی تحقیق با جنبه‌های کاربردی متعددی همراه گردد. 


مواد و روش‌ها 

منطقه موردمطالعه 

حوزه آبخیز گدارچای سومین حوضه Sop‏ حوضه آبریز دریاچه ارومیه در کشور ایران» بعد از زرینه‌رود و 
سیمینه‌رود است که به عرض‌های شمالی در موقعیت 1۵*۳1 تا "۳۷۹۱۰ و طول‌های شرقی در موقعیت 
"۶ تا "1۵0*۶۱ محدود است. وسعت این حوضه بالغ بر ۲۰۹۰ کیلومترمربع می‌باشد. رودخانه گدارچای 
پیه جیک از کوه‌های اشنویه سرچشمه و از منطقه جاشیران وارد دشت سولدوز شده و بعد از تغذیه آب‌های 
زیرزمینی و col‏ زمین‌های کشاورزی hle‏ طی طریق ۱۰۰ کیلومتر در روستای ممیند به dork yo‏ ارومیه 
منتهی می‌شود. رودخانه گدارچای دارای ۱۲ ایستگاه هیدرومتری است که مهم‌ترین آن‌ها ایستگاه‌های پی 
قلعه» پل Se ty‏ و پل بهراملو (جاده سنتو) است. بر روی رودخانه گدارچای سدی هنوز احداث نشده 
است؛ ولی سد حسنلو با حجم مخزن ٩۱۰‏ میلیون مترمکعب. که خارج از بستر این رودخانه است. از رودخانه 
گدارچای آبگیری می کند؛ و از لحاظ ویژگی‌های فیزیوگرافیکی حوضه گدارچای. دارای ۱۱ واحد کاری پل 
بهراملو سنتی محمد شاه سفلی» نقده بالچی, با یزید ay Obl‏ جیک. چپرآباده صوفیان, اشنویه. پی قلعه و 
دورود می‌باشد. پل بهراملو سنتو با مساحت ۲۰۲۹/۸۶ کیلومترمربع بیشترین وسعت و صوفیان با مساحت 
۹ کیلومترمربع کمترین وسعت را به خود اختصاص داده است. این حوضه در محدوده ارتفاعی ۱۳۹۳/۸۳ تا 
۲ متر از سطح Loe‏ قرار گرفته است. بیش‌ترین شیب hog‏ وزنی با مقدار 1۵/۲۲ Xe;‏ مربوط 4 
واحدکاری پی قلعه و کمترین با مقدار ۱۳/۰۱ درصد مربوط به Joly‏ کاری صوفیان می‌باشد (Yamani, & Abbasi,‏ 


)2020 در شکل )1( موقعیت جغرافیایی حوزه آبخیز گدارچای در ایران و استان آذربایجان غربی ارائه شده است. 


حوضه آبریز دریاچه ارومیه Eg‏ 
حوزه آبخیزگدارچای CIL‏ 


شکل ۱- موقعیت جغرافیایی حوزه آبخیز گدارچای- دریاچه ارومیه 
Fig.1. Geographical Location of the Godarchay Watershed - Lake Urmia‏ 


روش تحقیق 
در تحقیق حاضر از تصاویر ماهواره‌ای Sentinel2‏ با تفکیک مکانی E‏ متری و دارا بودن محدوده 
slab‏ متنوع (۱۳ ul‏ برای (Sentinel2‏ استفاده گردید. این تصاویر از اواخر سال ۲۰۱۵ میلادی از 


Copernicus cul.‏ آژانس فضایی اتحادیه اروپا دریافت گردید. اطلاعات تصاویر مورداستفاده در این 
مطالعه در حدول 0 ارائه شده است. 


حدول ۱- اطلاعات تصاویر ماهواره‌ای مورداستفاده در مطالعه تغییر و پیش‌بینی کاربری 


Table1- Information on Satellite Images Used in the Study of Land Use Change and 
Prediction 


قدرت تفکیک مکانی 


نام سنجنده تاریخ کسب داده 
Sentinel2‏ ۲۰۱۳/۰۷/۱۰ معادل ۱۳۹۵/۶/۲۲ چند طیفی ۲۰,۱۰ و ٩۰‏ متر 
Sentinel2‏ — ۲۰۲۰/۰۸/۲۰ معادل ۱۳۹۹/۰۶/۲ چند طیفی ۲۰,۱۰ و ٩۰‏ متر 
Sentinel2‏ — ۲۰۲۲/۰۹/۲۰ معادل ۱۶۰۱/۰۳/۳۰ چند طیفی ۲۰۱۰ و ٩۰‏ متر 


پیش پردازش تصاویر ماهواره‌ای 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف ۲۲۲۳ 


تصاویر E‏ ارسال لایه‌های AM‏ شده در dac‏ نرم‌افزارهای (ERDAS Imagine & ENVI)‏ انجام شد. 


-Y-Y-Y‏ در مرحله پردازش برای ایجاد فضای کار نرم‌افزاری» ورود داده‌ها و تصاوین ایجاد پروژه‌های 
موردنیاز. بعد از فراخوانی تصویر و ایجاد ترکیبی CHL‏ مناسب. که در این تحفیق از باندهای ۲۰۳۶۸ و ۱۱ 
استفاده شد. بعد از ایجاد ترکیب باندی سگمنت‌سازی (قطعه‌بندی) انجام شد. سگمنت‌سازی به مفهوم 
جداسازی پیکسل‌های تصویری در باندهای مختلف تصاویر ماهواره‌ای و یکپارچه‌سازی مجدد آن بر مبنای 
پارامترهای مختلف نظیر شباهت‌های بازتابشی» هندسی» تراکمی» فشردگی, ارتباط ذاتی آن‌ها باهمدیگر از 
Gb‏ اعمال الگوریتم‌های تقسیم‌بندی متعدد تعریف می‌گردد (رسولی۱۳۸۷). نرم‌افزار eCognition‏ 
Developer 9.5‏ جهت طبقه‌بندی راه‌های مختلفی در پیش روی کاربر قرار می‌دهد. طبقه‌بندی‌هایی که We‏ 
در این زمینه و با استفاده از این نرم‌افزار انجام گرفته به‌صورت مراحل سگمنت‌سازی» نمونه‌برداری و 
طبقه‌بندی است. ازآنجایی که در مرحله طبقه‌بندی شیء گرا تعریف شیءها یا سگمنت‌ها است» به همین منظور 
برای سگمنت‌سازی, مقدار OT‏ ۰/۳ و ۰/۷ به ترتیب برای پارامترهای مقیاس. شکل و فشردگی مطلوب بودن 
در نظر گرفته شد. جهت تفکیک pe‏ پوشش گیاهی از شاخحص (NDWI')‏ و (NDVI?)‏ استفاده شد. رابطه 
شاخص‌ها برای ماهواره Sentinel2‏ در abl,‏ ۱ و Y‏ بیان شده است (Rasuly Pirouzian, Chung,‏ 
mE Moharrami, & Derafshi, 2015)‏ 


NDVI Sentinel2 = (Near Infraredb8 — REDb4 / Near Infraredb8 + REDb4) (1) 


NDWI Sentinel2 = (Greenb3 — Near Infraredb8 / Greenb3 + Near Infraredb8) (2) 

پس از ایجاد شاخص‌هاء با توجه به هدف تحقیق که آشکارسازی تغییرات کاربری اراضی است» کلاس‌های 
کاربری در نظر گرفته‌شده شامل کاربری پهنه آب. مسکونی مرتع» نمک حاشیه دریاچه ارومیه. باغات و 
اراضی آبی» اراضی دیم و خاک است. در طبقه‌بندی تصوير با روش شیءگرا از الگوریتم نزدیک‌ترین همسایه 
NNC‏ استفاده شد و متناسب با خصوصیات بصری و به کمک شاخص‌ها نمونه برای هرکدام از کلاس‌ها 
انتخاب و lulu‏ طبقه‌بندی برای هر یک از کلاس‌های کاربری اراضی تعریف شد و نقشه‌های کاربری 
اراضی برای Y‏ دوره با روش شیءگرا استخراج گردید و مساحت کاربری‌ها طبقه‌بندی شده محاسبه گردید. 

در مرحله بعد پیش‌بینی تغییرات کاربری به کمک مارکوف و )ما رکوف انجام شد. 


1. Normalized Difference Water Index 
2. Normalized Difference Vegetation Index 


۱۶۰۳ تابستان‎ Y شماره‎ Wile damus نشریه جغرافیا و مخاطرات‎ ٤ 


با توجه به توانایی الگوریتم تلفیقی زنجیره مارکوف و سلول‌های خودکار از این مدل برای پوشش سطح 
زمین و پیش‌بینی کاربری در آینده استفاده می‌شود. این مدل با تحلیل یک جفت از نقشه‌های کاربری زمین با 
استفاده از زنجیره مارکوف» ماتریس احتمال انتقال ماتریس مساحت انتقال یافته" و چند تصویر احتمال 
شرطی ایجاد می‌گردد (بسته به تعداد کاربری‌ها)» ماتریس احتمال انتقال, احتمال انتقال هرکدام از کاربری‌ها 
به کاربری So‏ را نشان می‌دهد )2018 (Mirakhorlo, & Rahimzadegan,‏ « این مدل پیش‌بینی آینده را 
با ایجاد یک ماتریس انتقال احتمال تغییرات کاربری بین سال اول (To)‏ و دوم (T1)‏ انجام می‌دهد. محاسبه 


پیش‌بینی تغییرات کاربری زمین از طریق رابطه (۳) به دست می‌آید: 


S =(T, +T,)=Pi *S (t) (۳ 
(Hamad & Kolo, 2018) احتمال انتقال می‌باشد‎ 

اساس تولید ماتریس احتمال انتقال این است که می‌توان از شرایط قبلی کاربری Os‏ برای پیش‌بینی استفاده 
کرد برای One‏ ماتریس احتمال انتقال از رابطه )£( استفاده می‌شود: 

Pu ۰ ۳۱2۰۰۰ (£)‏ 
Py Pon Pas‏ 
که در Pij ol‏ احتمال انتقال را از اولین سال G)‏ تا دومین سال G)‏ نشان می‌دهد & (Mirakhorlo,‏ 


Pij = (0x Pij 21) 


Rahimzadegan, 2018) 

مدل‌سازی تغییر کاربری اراضی ۲۰۲۲ با مدل CA‏ مارکوف 

جهت تولید نقشه پیش‌بینی کاربری اراضی YY‏ ابتدا به‌وسیله مدل مارکوف از نقشه کاربری اراضی ۲۰۱۱ و 
۶۰ برای agi‏ ماتریس احتمال انتقال و ماتریس مساحت انتقال جهت صحت‌سنجی نقشه پیش‌بینی استفاده 
شا وروی فان مان کف ویر دی ا ce Ly a‏ پم ای کی وی مال 
cil aiat‏ بش ۲ سال) می اک 

مرحله پس‌پردازش ورود لایه‌های رقومی» فرمت‌سازی تصاویر طبقه‌بندی شده» ارزش‌گذاری کلاس‌های 
کاربری اراضی» تولید لایه‌های پیش‌بینی کاربری اراضی» ارزیابی نتایج نهایی و ارسال نتایج به محیط نرم‌افزار 
5 ترکیب و تحلیل هدفمند لایه‌ها انجام می‌شود و نتایج نهایی استخراج گردید. در شکل (Y)‏ 


فرآیند کلی انجام تحقیق ارائه شده است. 


دریافت و مدیریت تصاویر ماهواره‌ای 


Sentinel-2 


1. Transition Probabilities 
2. Transition Areas 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف YYO‏ 


پیش - پردازش تصاویر 
(ERDAS Imagine & ENVI)‏ 


انتقال تصاویر به محیط نرم افزار 


eCognition Developer 9.5 


تولید شاخص‌های طیفی اعمال شروط. آستانه گذاری و قطعه‌بندی جند مقیاسی- تفکیکی 
فیلترهای قانونمند تطبیقی و شئ گرا 
TUE din d NDVI and NDWI‏ 
طبقه‌بندی نظارت‌شده به روش NNC‏ محیط نرم‌افزار 


تولید نقشه‌های کاربری اراضی ۲۰۱۹ ۲۰۲۰ و ۲۰۲۲ 


EDRSI -TerrSet ۲۰۲۲ کاربری اراضی‎ 


= 


بررسی روند تغییرات 
پیش‌بینی کاربری اراضی ۲۰۲۲ 
حوادث IS‏ شته 


انتقال یافته‌های تحقیق به محیط GIS‏ 


(تحلیل‌های پس-پردازش و کارتوگرافی) 


بررسی جنبه‌های 


کاربردی مدل‌های حاصله 


شکل ۲- نمودار جریانی مراحل انجام تحقیق 
Fig.2. Flowchart of research steps‏ 


a i ۳۳1‏ جغرافیا و مخاطرات محیطی. جلد۱۳. شماره ۲ تابستان ۱۶۰۳ 


نتایج و بحث 


پردازش تصاویر به روش‌های شیء گرا پیکسل‌های تصویر ماهواره‌ای (تشکیل‌دهنده باندهای تصویری) با 
اعمال الگوریتم‌های سگمنت‌سازی (Segmentation)‏ و با توجه به میزان بازتابش» شکل» تراکم و phe‏ 
ارزش‌های زمین-آماری به شیءهای تصویری تبدیل گردید با اعمال روش‌های ایندکس‌سازی. طبقه‌بندی 
تصاویر انجام گردید. در این تحقیق از ترکیب باندهای ۲ A . Y‏ و ۱۱ تصویر 560110612 استفاده شده 


است. نتایج ترکیب باندی تصویر» سگمنت‌سازی و طبقه‌بندی تصویر در شکل (۳) ارائه شده است. 


Image band composition 
تر کیب باندی تصویر‎ 


Segmentation 
سگمنت سازی (قطعه بندی تصویر)‎ 


Image classification 
طبقه‌بندی تصویر‎ 


شکل ۳- تر کیب باندی. سگمنت سازی و طبقه‌بندی تصویر 


Fig.3. Band composition, segmentation and image classification 


بعد از طبقه‌بندی تصاویر به روش NNC‏ در هفت کاربری پهنه‌آبی. CEL SE‏ و اراضی‌آبی» اراضی دیم نمک 
حاشیه دریاچه ارومیه» مرتع و مسکونی نقشه‌های کاربری اراضی سال ۲۰۱۰ با دقت ٩۱‏ درصد. سال ۲۰۲۰ با دقت 
۵ درصد و سال ۲۰۲۲ با دقت ٩۳‏ درصد تهیه گردید. مساحت هر کاربری محاسبه شد. شکل )£( نقشه‌های کاربری 
اراضی سال Ye Y‏ ۲۰۲۰ و ۲۰۲۲ نمایش می‌دهند. طبق نتایج برای سال ۲۰۱۳ مساحت پهنه SE ul‏ باغات و 
"PP‏ اراضی دیم. نمک حاشیه دریاچه ارومیه. مرتع و مسکونی به ترتیب YY AAV OVA VAY‏ ۳۵۱/۱۹ 
۵ ۰/۷۰ و ۲۲/۵۲ کیلومترمربع به دست آمد. برای سال ۲۰۲۰ مساحت پهنه‌آبی» خاک» باغات و اراضی آبی 


اراضی Kas qu‏ حاشبه دریاجه ارومیه. مرتع E‏ مسکونی به ترئیب ۹ ۵۳۵/۵ ۵۰۱۸/۳۶ (Y/YY OVAN‏ ۳۵/۶۵" 


سال سيزدهم T‏ و همکاران. پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف YYV‏ 


و ۲۶/۰۵ کیلومترمربع» همچنین در سال ۲۰۲۲ مساحت برای کاربری‌های پهنه‌آبی. خاک CEL‏ و اراضی cel‏ اراضی 


کیلومترمربع بود. 
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شکل -٤‏ نقشه کاربری اراضی برای سال‌های موردبررسی در حوزه آبخیز گدارچای-دریاچه ارومیه 
Fig.4. Land Use Maps for the Years Under Study in the Godarchay Watershed - Lake‏ 


Urmia 

صحت‌سنجی مدل CA‏ مار کوف 
به‌منظور صحت‌سنجی نقشه طبقه‌بندی ۲۰۲۲ با نقشه پیش‌بینی ۲۰۲۲ مقایسه گردید. مقادیر Kno‏ برابر با 
٩‏ درصد plp Klocation‏ با ۰/۹۲ درصد. KlocationStrata‏ برابر با ۰/۹۲ درصد و کاپای استاندارد 
(Kstandard)‏ برابر با ۰/۸۲۳ درصد محاسبه گردید و با توجه به دقت بالای مدل در پیش‌بینی نقشه تغییرات 
سال ۲۰۲۲ اقدام به استخراج نقشه تغییرات ۲۰۲۸ شد. در این راستاء در مطالعه‌ای )2017 (HUA,‏ از تصویر 
طبقه‌بندی‌شده و تصویر پیش‌بینی‌شده به‌وسیله مدل CA‏ مارکوف برای صحت‌سنجی استفاده شده است 
ضریب کایای استاندارد برابر با ۰/۸۶۳ کایای موقعیت برابر با ۰/۸۱۳ و کاپای موقعیت استانذار برابر با 


Salada] به‌دست‎ ۳ 


۱۶۰۳ تابستان‎ Y شماره‎ Whe نشریه جغرافیا و مخاطرات محیطی‎ YYA 


ماتریس مساحت انتقال و ماتریس احتمال انتقال جدول (۲) و (۳) با استفاده از نقشه طبقه‌بندی سال VAT‏ و 
۰ به‌دست آمد. ماتریس مساحت انتقال ۲۰۱-۲۰۲۰ به‌عنوان ورودی به مدل CA‏ معرفی شد تا نقشه 


پیش‌بینی سال ۲۶۲۲ به‌دست LI‏ 


حدول ۲- ماتریس مساحت انتقال یافته تغییر کاربری (۲۰۱۹-۲۰۲۰) با استفاده مدل زنحیره مار كوف بر حسب 


کیلومترمربع 


Table2- Land Use Change Transition Area Matrix (2016-2020) Using the Markov Chain 
Model in Square Kilometers 


کاربری اراضی پهنه آبی خاک مناطق مسکونی باغات و اراضی آبی مرتع = نمک حاشیه دریاچه 

JM ۰/۳۹ ۱/۹۹ YN . AN 1/44 پهنه آبی‎ 
Y ۱۹۳/۵ ۲ ADE NW OAV/Y \ YA خاک‎ 

۰/۰١ ۲/۷ Y/AÉ ۰/۸۹ vv ۹۳ uy مناطق مسکونی‎ 
11۹۵ Yav 33711 010/07 ۱/۳۳ Vo/VA VAL و اراضی آبی‎ SEL 
۳/۹۵ ۱۰۱/۷۳ ۵ YAV/00 ۱/۳۳ ۳۲۳/۳۹ vor مرتع‎ 
aes V/A 2/17 0/4۸ ۱/۳۹۵ ۱۳۳/۰۳ ۱/۳۵ ديم‎ 
uw . T ۰/۱ . ۱/۱ a نمک حاشیه دریاچه‎ 


براساس نتایج حاصل از جدول(۲)» 1۰/۹۹ کیلومترمربع از پهنه‌آبی در سال ۲۰۱۳ در سال ۲۰۲۰ در همان 


باغات و اراضی آبی و ۲/۹۱ کیلومترمربع به کاربری مرتع در سال ۰ تبدیل خواهد شد. 


حدول۳- ماتریس Jul Jul‏ (پیش بینی شده برای YY‏ با استفاده از طبقه بندی ۰ (Y YM- Y‏ بر حسب درصد 


Table3- Transition probability matrix (predicted for 2022 using 2016-2020 classification) 
in percentage 


کاربری اراضی پهنه آبی خاک مناطق مسکونی باغات و اراضی آبی . مرتع = نمک حاشیه دریاچه 
پهنه آبی ۷۱/۳۸ Ji ۳/۶ ۲ . YA‏ ۰/۸ 
خاک Zu VY ۳/۵ “/V £aY VAY‏ 

مناطق مسکونی 13 1/۳ 90/۲ ۳ ۱۳/۱ vit‏ 

UA 1/۸ ۲/۱ ۳/۷ ZA ۹/۲ Nt و اراضی آبی‎ SEL 
A AA Y4/4 YY (A ۰/١ YvY ۰/۵ مرتع‎ 
۰/۱ WY YYA a/y Y ۲۳/۸۷ ۹ ديم‎ 
۰/۱ . VE ۰/۱ ` NY ۷۹/۶۲ نمک حاشیه دریاچه‎ 


TRA‏ خواهد ماند. هم ARE edo yj ۲۰/7 ure‏ درصد. Y‏ درصد. ۶ درصد و ۱/۶ درصد از کاربری 


پهنه‌آبی در سال ۲۰۱ به کاربری‌های خاک مناطق مسکونی CEL‏ و اراضی coal‏ مرتع و ديم تبدیل خواهد 


شال atis cud‏ ایرانی و همکاران. پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف ۳۳۹ 


شد. در این ماتریس ردیف‌ها نشان‌دهنده سال شروع (VU‏ و ستون‌ها نشان‌دهنده سال پایان (۲۰۲۲) 
مدل معرفی مت علاوه بر ماتریس مساحت انتقال نقشه طبقه‌بندی کاربری اراضی 1 به‌عنوان تصویر eb‏ 
و تعداد دوره چهار سال به‌عنوان ورودی به مدل تعریف شد و نقشه پیش‌بینی ۲۳ شکل )0( تهیه شد. 
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شکل ۵- نقشه پیش‌بینی کاربری اراضی سال ۲۰۲۲ 

Fig.5. Land Use Prediction Map for 2022‏ 
بعد از دستیابی به دقت مورد نیاز برای نقشه پیش‌بینی سال ۲۰۲۲ پیش‌بینی برای سال ۲۰۲۸ انجام شد. برای 
ضریب Lis‏ نشان‌دهنده Ol x‏ دقت پیش‌بینی مدل می‌باشد مقدار این ضریب بین و ۱ می‌باشد که در این 
تحقیق دقت CYL‏ ۰/۸ به‌دست آمد؛ که هرچه‌قدر مقدار OF‏ نزدیک به یک باشد نشان‌دهنده دقت بالای مدل 
می‌باشد. نتایج این تحقیق با نتایج )2017( Yirsaw, Wu, Shi, Temesgeh, & Bekele,‏ در مطالعه‌ای ضریب 


در آینده می‌باشد مطابقت دارد. برای به‌دست آوردن ماتریس مساحت انتقال ابتدا به‌وسیله مدل مارکوف 


.۳۳ 4 جغرافیا و مخاطرات محیطی. جلد۱۳. شماره Y‏ تابستان ۱۰۳ 


TAE‏ مال WV‏ و sles 0١‏ قد کا مار س ات اتفال seal‏ کد و ادان زر ود به 
مدل CA‏ مارکوف برای پیش‌بینی سال ۲۰۲۸ معرفی شود. فاصله زمانی بین دو تصویر 1 سال و فاصله زمانی 
ely‏ میتی o‏ سال Vas‏ سال oe‏ کته شف ر ail‏ سفق ماس Sery (EH els‏ 
۵ نقشه پایه (۲۰۲۲) وتعداد کل تکرارها به‌عنوان ورودی‌های مدل CA‏ مارکوف بود. تعداد کل تکرارها 
بر اساس زمان تعیین‌شده برای پیش‌بینی تغییرات کاربری اراضی می‌باشد (Sarabuddin Mondal, Sharma,‏ 
Kappas, & Garg, 2019)‏ که در این مطالعه تعداد دوره‌ها 1 سال انتخاب شد. ماتریس مساحت انتقال و 
مار اف اال Ud‏ سیون OY‏ تفت ud‏ از یه d‏ سا VO‏ ۶ ۲۸۲۲ سس 


al 


حدول £~- ماتریس مساحت انتفال L (YA Y YY) ash‏ استفاده مدل زنحیره‌ای مار كوف بر حسب کیلو مترمربع 
Table4- Transition Area Matrix (2016-2022) Using the Markov Chain Model in Square‏ 


Kilometers 
نمک حاشیه دریاچه ارومیه‎ E e^ کاربری اراضی آب خاک مسکونی باغات و اراضی آبی‎ 
NO ۰/۳۸ (A ۱۳:۳ TIE AE ۳۳/۷۲ پهنه آبی‎ 
۱۳/۳۷ ۹7/۷۱ 14/۰0 7۱/10 wav Mv 0/۷ خاک‎ 
v/s) ۳/۳۹ ۱۳ Y/*^ ۳۰/۸ VAN ۵ مسکونی‎ 
£/AN oY/í- YV/YN VY+ /0V ۱/۹۰ ۳۹/۳۰ 0/0۷ al باغات و اراضی‎ 
£/40 ۹/7۱ ۹/٤۱ YAVA \/VA 1۸/10 ۳/۲ مرتع‎ 
۳۷ YY\/00 ۱۳۳/۳۰۶ ۱۳3۹/۷۵ EV ۱۷ ۳۹ E 
۱/۳۸ VE ۰ VES ۰/١ ۱۳/۰ ۰۵ نمک حاشیه دریاچه ارومیه‎ 


براساس gold‏ حاصل از جدول )6 WAY‏ کیلومترمریم از کاربری آب در سال ۲۶۱۸ در سال ۲۰۲۲ دز 
همان کاربری آب BL‏ خواهد ماند. هم‌چنین. ۱۰/۲۰ کیلومترمربع از کاربری آب به خاک» ۰/۰۶ کیلومترمربع 
به کاربری مسکونی» VET‏ کیلومترمربع به کاربری باغات و اراضی آبی:۰/۷۱ کیلومترمربع از کاربری آب به 
مرتع» ۰/۲۸ کیلومترمربع به کاربری ديم و ۰/۷۵ کیلومترمربم به نمک حاشیه دریاچه ارومیه در سال ۲۰۲۲ 
des‏ خواهد شد. ماتریس احتمال انتقال. احتمال اینکه هر نوع کاربری اراضی را در یک دوره زمانی معین 


در آینده مربوط به هر پیکسل در تصوير باشد گزارش می‌دهد. 


حدول 0- ماتریس احتمال انتقال (پیش بینی شده برای L Y* YA‏ استفاده از طبقه‌بندی (Y*YM- Y YN‏ بر حسب 


درصد 


Table5- Transition probability matrix (predicted for 2028 using 2016-2022 classification) 
in percentage 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف YYY‏ 


کاربری اراضی آب خاک مسکونی ‏ باغات و اراضی آبی مرتع = نمک حاشیه دریاچه ارومیه 

۲ ۷۵ IE ۳/۸۵ ۱ YV/££ W PE 

۳۳ \V/0 ۵۳ ۸ vy) VYY ové: ۰/1: خاک‎ 

۳ VAA MUNI ۳ "۰۰۲۵ ۰ ۰ ۰ ۰ ۰ ۲ مسکونی‎ 

باغات و اراضی آبی ۳ Y/£^‏ ۰/۸ ۳/۷ ۱۹/۹۰ ۵ 23 

۰/۹ ۹/4۲ AA) YV/VO ۰/۸ A10 AS مرتع‎ 

TR YV/ÁÉ YY ۱۹/۳ ۰:۲ UIT ۸ ديم‎ 

۸/۳۹ ۷ v/s) WAR! IN AV/04 ۳/۵۵ نمک حاشیه دریاچه ارومیه‎ 


با توجه به نتایج جدول )0( WA‏ درصد از کاربری پهنه آبی در سال ۲۰۱۳ در همان کاربری پهنه آبی در 
سال ۲۰۲۸ باقی خواهد ماند. هم‌چنین» ۲۷/4۶ درصد. ۰/۱۱ درصد. ۳/۸۵ درصد. ۲/۰۶ درصد. ۰/۷۵ درصد 
و Y‏ درصد از کاربری پهنه‌آبی در سال ۲۰۱۳ به کاربری‌های اراضی SE‏ مسکونی. باغات و اراضی coal‏ 
مرتع» کشاورزی ديم و نمک حاشیه ambos‏ ارومیه تبدیل خواهد شد. در این ماتریس ردیف‌ها نشان‌دهنده 
سال شروع (VND)‏ و ستون‌ها نشان‌دهنده سال پایان (۲۰۲۸) می‌باشد. خروجی مدل مارکوف (ماتریس 
مساحت انتقال) به‌عنوان یکی از ورودی‌های مدل CA‏ مارکوف به مدل معرفی شد. علاوه بر ماتریس 
مساحت انتقال نقشه طبقه‌بندی کاربری اراضی ۲۰۲۲ به‌عنوان تصویر پایه» تعداد دوره 7 سال به‌عنوان ورودی 
به مدل تعریف شد و نقشه پیش‌بینی ۲۰۲۸ تهیه شد. در تهیه نقشه پیش‌بینی سال Y* YA‏ از ماتریس احتمال 
انتقال و ماتریس مساحت انتقال به‌دست آمده سال‌های ۲۰۱-۲۰۲۲ استفاده شد. در بازه زمانی سال‌های 
۲۰۲۲-۹۲ بیش‌ترین مساحت انتقال در کاربری‌های ذکرشده به ترتیب مربوط به SEL SE‏ و اراضی 
آبی دیم مرتع» مسکونی, پهنه‌آبی و نمک حاشیه دریاچه ارومیه بود. ماتریس مساحت انتقال» ميزان تغییر در 
تعداد پیکسل‌ها را از یک طبقه کاربری به طبقه دیگر در دوره مشخص را نشان می‌دهد. دلیل بالا بودن 
مساحت انتقال در کاربری مرتع و خاک را می‌توان افزایش اراضی آپی و باغات دانست. نتایج OUS‏ می‌دهد 
که مساحت کاربری‌های کشاورزی (کشاورزی آبی و باغات) به سرعت افزایش پیدا کرده است به‌طوری‌که 
مساحت کاربری کشاورزی آبی از ۲۰/۷۸ درصد کل مساحت حوضه در سال ۲۰۱۲ به بیش ۳۶/۵۷ درصد در 
سال ۲۰۲۸ افزایش پیدا خواهد کرد؛ و همچنین مساحت کاربری مرتع و خاک به ترتیب ۲۸/۹۲ و ۳۱/۵۷ 
درصد کل مساحت حوضه در سال ۲۰۱۲ به ۲۳/۵۷ و ۲۱/۳۰ درصد در سال ۲۰۲۸ کاهش پیدا خواهد کرد 
که با نتایج )2019 (Kangabam, Selvaraj, & Govindaraju,‏ در ارزیابی تغییرات کاربری اراضی و پوشش 
گیاهی برای دلتای نیل مصر در طول YA‏ سال» بیان نمودند گسترش مداوم اراضی کشاورزی و مناطق شهری» 
هم‌زمان با کاهش در پهنه آب و مراتع رخ‌داده است مطایقت دارد. نقشه کاربری پیش‌بینی شده سال ۲۰۲۸ به 
روش CA‏ مارکوف با ضریب LIS‏ ۰/۸۷ تهیه گردید؛ که در شکل OU)‏ آمده است. 


۳۳۲ نشریه جغرافیا و مخاطرات dons‏ ¢ جلد۱۳» شماره «Y‏ تابستان VEY‏ 
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Y* YA شکل - نقشه پیش بینی کاربری اراضی سال‎ 
Fig.6. Predicted Land Use Map for 2028 


آشکارسازی تغییرات برای درک تغییرات منطقه» گسترش و کاهش مساحت کلاس‌ها در سال‌های ۲۰۱7 


۲ و۲۰۲۸ به کیلومترمربع و درصد در جدول (1) ارائه شده است. 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف ۲۲۲ 


جدول *- مساحت کاربری اراضی به کیلومترمربع و درصد 
Table6- Land use area in square kilometers and percentage‏ 


Y YA ۲۰۳۲ PY ۲۰۱۹ سال‎ 

کلاس کاربری اراضی کیلومترمربع درصد کیلومترمربع درصد کیلومترمربع درصد کیلومترمربع درصد 
پهنه آبی VY YV/Y' +/40 MAY‏ ۳/۳ ۱/۳۲ ۳/۱ ۰/۹۹ 

Wir ££0/0 YE/va ۰9۸/۵ Yo/M 00/0 1/0۷ 11/۰۹ خاک‎ 

\/£) ۳۹/۶ VE ۳۸/۲ ۱/۵ ۱/۳۵ ۱/۳۸ YY/oY مسکونی‎ 
Y£/oV VYY/A ۳۹/۹۶ ۳3/۲ Y£/Y OAS YUM £YY/o- باغات و اراضی آبی‎ 
YY/oV 1۹۳/۸ YE/AY AK ۳۰/۳۹ 10/6۹ YA/AY ۷۰ مراتع‎ 

۱/۹۳ Yvo ۱۳/۹۳ YV£/A ۱۷/۲ 0/41 ۱3/۸۰ ۱۳9۵۹ اراضی دیم‎ 

VAR! o/t AY Y/^ JM ۲/۳۱ ۲ ۵ نمک حاشیه دریاچه ارومیه‎ 


با توجه به جدول O)‏ بیشترین تغییرات از کل مساحت منطقه. مربوط به کلاس اراضی CEL‏ و اراضی آبی 
(کشاورزی آبی) می‌باشد به گونه‌ای که مساحت این کلاس از ۶۳۲/۵۰ کیلومترمربع در سال ۲۰۱۹ به 1۲۷/۲ 
کیلومترمربع در سال ۲۰۲۲ رسیده است به بیان دقیق‌تر مساحت دو کاربری خاک و مرتع مقدار کاهشی و 
اراضی دیم. مسکونیء باغات و اراضی آبی و نمک حاشیه دریاچه ارومیه مقدار افزایشی داشته است؛ و 
پیش‌بینی می‌شود در سال ۲۰۲۸ مساحت کاربری خاک و مرتع کاهشیء و مساحت کاربری مسکونی, باغات 
و اراضی qul‏ ديم و نمک حاشیه دریاچه ارومیه افزایش LL‏ درنهایت کاربری‌های ديم مسکونی و نمک 
حاشیه دریاچه ارومیه با مساحت NVA‏ ۱/۰۸ و ۰/۰۰۲ درصد کل مساحت در سال Y» YA‏ به VAY‏ ۱/۶۱ و 
٩‏ درصد کل مساحت در سال ۲۰۲۸ خواهد رسید. افزایش مساحت کاربری‌های کشاورزی در حوزه 
آبخیز گدارچای بخصوص دو کاربری CEL‏ و زمین‌های کشاورزی آبی. باعث افزایش مصرف منابع آب 
خواهد شد. دلیل بالا بودن احتمال انتقال کاربری اراضی آبی و CEL‏ را می‌توان گرایش مردم به سمت ایجاد 
باغ و درآمدهای حاصل از باغ‌داری در منطقه موردمطالعه دانست. براساس پیش‌بینی مدل CA‏ مارکوف سطح 
کشاورزی آبی و باغات در مجاور منابع آبی و نیز در مجاور اراضی کشاورزی آبی موجود افزایش خواهد 
یافت که دلیل آن را می‌توان در با توسعه زراعت آبی و استفاده از روش‌های آبیاری قطره‌ای و غیره توسط 
کشاورزان مرتبط دانست. علاوه بر این» دلیل افزایش مساحت باغات را می‌توان علاقه مردم به داشتن باغ 
شخصی و درآمد ناشی از محصولات باغی دانست. نتایج این تحقیق با نتایج Aburas, Abdullah, Ramli,‏ 
Ash'aari, & Ahamad, (2018)‏ اذعان کردند که ناحیه‌های کشاورزی از ناحیه‌های شهری در ۱۰ سال 
آینده تأثیر خواهد گرفت. که این امر ناشی از توسعه اقتصادی خواهد oy‏ مطابقت دارد. همچنین نتایج این 
Ges‏ با نتایج )2021( Rasouli, ۸۵8270۷۵, & Safarov,‏ در مطالعه‌ای شناسایی و ارزیابی تغییرات 
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LLU/LC‏ استفاده از تصاویر ماهواره‌ای Sentinel2‏ سال ۲۰۱۹ و ۲۰۲۱ با بکار گیری تکنیک‌های پیشرفته 
شیءگرا انجام گردید. درنهایت» برای پیش‌بینی انواع SLU‏ مدل زنجیره مارکوف خودکار سلولی CA‏ 
استفاده شد. نتایج نشان داد تغییرات کاهشی در انواع LU‏ عمدتاً در جنگل‌ها و مراتع» و روند افزایشی قابل 
توجهی در زمین‌های ub‏ و رها شده وجود خواهد داشت. دقت کلی ۰/۹۵ و ضریب LE‏ نزدیک به ٩۶‏ 
درصد بود که تأیید کرد تغییرات bE‏ توجهی بر نقشه‌های LU‏ منطقه موردمطالعه وجود داشته است مطابقت 


دارد. روند تغییرات کاربری اراضی در شکل (V)‏ آمده است. 
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شکل ۷- روند تغییرات کاربری اراضی به روش شیءگرا 
Trends in Land Use Changes Using the Object-Oriented Method‏ ۲۱۵۰7۰ 
کاربری‌ها (تبدیل کاربری مرتع و خاک به کل اراضی آبی و EL‏ دیم و مسکونی) را نشان می‌دهد که با 


نتایج )2019 (Birhanu, Masih, van der Zaag, Nyssen, & Cai,‏ در تحقیق بر روی حوضه PHI‏ گومارا! اتیویی 


2. Gomara 


سال سیزدهم ایرانی و همکاران پیش‌بینی تغییرات کاربری اراضی با استفاده از روش‌های پردازش شیءگرا و زنجیره مارکوف YYO‏ 


نشان دادند در مدت ۳۹ سال TUE.‏ مساحت مراتع als‏ یافته و در مقابل مساحت اراضی کشت شده 


افزایش یافته و همین امر باعث کاهش حجم آب‌های جاری حوضه آبریز شده است مطابقت دارد. 


نتیجه گیری 


در این تحقیق برای آشکارسازی و تحلیل تغییرات سال ۲۰۱۳ و ۲۰۲۸ از تصاویر سال‌های ۰۲۰۱7 ۲۰۲۰ و 
esl ۷۲‏ لب رای ll‏ مدل 6 نیسای 238 edil ce Y Lo‏ عاص 3l‏ 
طبقه‌بندی تصویر ماهواره‌ای همان سال مقایسه شد. ضریب SUIS‏ به‌دست آمده بیانگر قابلیت بالای مدل CA‏ 
مارکوف برای شبیه‌سازی تغیبرات کاربری اراضی در حوزه آبخیز گدارچای بود بعد از تهیه نقشه‌های کاربری 
سال ۲۰۱۹ ۲۰۲۰ و ۲۰۲۲ و دستیانی به دقت مورد to fe abe cele JU‏ مال ۲۰۲۳ با استفاده از de‏ 
CA‏ مارکوف. پیش‌بینی برای سال ۲۰۲۸ انجام شد و نقشه تغییرات کاربری اراضی برای سال ۲۰۲۸ تهیه شد. 
ضریب کاپای ۰/۸۷ به‌دست آمده بیانگر قابلیت بالای مدل CA‏ مارکوف برای پیش‌بینی تغییرات کاربری 
اراضی در حوزه آبخیز گدارچای بود. نتایج OLE‏ دهنده افزايش اراضی آبی و باغات» مسکونی و ديم در 
سال‌های آتی می‌باشد. براساس نتایج تحقیق, می‌توان در راستای کنترل تغییر کاربری اراضی متناسب با قابلیت 
اراضی و اتخاذ استراتژی‌های کاربردی در حفظ اراضی طبیعی اقدام نمود. قابل ذکر است که نتایج حاصل از 
تفسیر سناریوهای مختلف کاربری اراضی می‌تواند در پیش‌بینی اثرات تغییر کاربری اراضی بر مژلفه‌های تولید 
رسوب. سیل‌خیزی و نیز اثرات اکولوژیک مورد استفاده قرار گیرد. Ab‏ اشاره شود که کاربری اراضی به‌عنوان 
عامل موازنه مولفه‌های اکولوژی. هیدرولوژی و اقتصاد به‌شمار می‌رود که تغییرات آن می‌تواند در ارزیابی 


اثرات وضعیت اقتصادی و نیز معیشت آبخیزنشینان اسفاده شود. 
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